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[ We^ compute the pair annihilation cross section of light dark matter scalar particles 

into two photons, and discuss the detectability of the monochromatic line associated 
with these annihilations. 
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O 1 1 Introduction 

The need for Cold Dark Matter (DM) to describe and understand how structures in the 
universe hold together, has become increasingly pressing with the impressive amout of 
observational data collected in the last ten years. Alternatives to DM, like modifications 
of gravity, are being put to critical and maybe fatal test by recording maps of gravitational 
lensing. Indeed, the separation recently observed in colliding clusters between the (maybe 
modified) gravitational deviation of light and the normal matter that causes it, seems very 
contrived with modified gravity, and very natural if collisionless DM is the main source 
of gravity. Questions about the nature of DM and its non-gravitational detection are 
therefore more relevant than ever. 

In this context, the precise determination by INTEGRAL / SPI^ of the characteristics 
of the 511 keV line emitted in our galaxy is intriguing. Indeed, it implies wihtout any 
doubt that the central bulge of our galaxy is a strong source of positrons. Astronhysical 
sources like Low Mass X-ray Binaries (LMXB) and Type la Supernovae (SNIA)^ cannot 
naturally explain why this source is at the same time steady, extended and absent in 
the disk. On the contrary, production of positrons through DM annihilation is naturally 
steady and concentrated in the bulge where the DM density increases: a fit of the needed 
DM density profile can even be attempted^, yielding a reasonable NFW profile p(r) ~ r _1 
for DM annihilation at rest, in opposition to a less reasonable p ~ r~ 2 for DM decay. 

In order to m aximize the electron-positron annihilation channel, such DM must be 
light (LDM)f^, at least below muon pair threshold: m dm < 100 MeV. More constrain- 
ing upper bounds can be obtained by a careful study of final state radiation processes 
(dm dm* — > e + e~ 7) and positron anihilation in flight, both producing continuous gamma 
ray spectra which increase with m dm . From the first, m dm < 35 MeV is obtained^, and 
m dm < 20 MeV from comparing the second with error bars on the measured spectra. On 
the other hand, m dm should be higher than 2 MeV to avoid spoiling nucleosynthesis^, 
and higher than 10 MeV if there is a significant coupling to neutrinos which can alter 
supernova explosions^, but such is not necessarily the case. 
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A fairly unique viable model satisfying all the above constraints contains scalar DM 
particles with m dm ~ 10 MeV, annihilating at rest in the galactic bulge into e + e~ pairs 
via t-channel exchange of heavy (> 100 GeV) fermions F e . Given the large local DM 
abundance inferred from the rotation curve, the annihilation cross-section yielding the 
observed positron source is however too small to explain a correct relic density inferred 
e.g. from cosmic microwave background measurements. A further light vector particle can 
then be invoked to mediate an s-channel annihilation process: being velocity dependent, 
this process becomes dominant in the early universe and can independently be adjusted 
to the relic density. 

If correct, such a model^ would profoundly alter the road to unification in particle 
physics. It therefore seems important to look for other experimental cross-checks. The 
simplest and most convincing one would be the discovery of another gamma ray line, from 
the process dm dm* — ► 77. In the following, we show^ this line is anavoidable in such a 
model, estimate its intensity, and discuss its observability. 



2 Dark matter annihilation cross section into two photons 

The model considered is specified by the Lagrangian C = ^F e ( c rPL + ciPn)ipe(fidm + h.c 
where Pr,l are the chiral projectors (1 ± 75) /2. The relevant annihilation diagrams are 
box-diagrams containing 1, 2 or 3 heavy fermions F e . Assuming that dm 7^ dm* (which 
fixes the circulation of arrows), there are 6 diagrams, taking into account permutation of 
the 2 photon external legs. 

From naive power counting, each box is logarithmically divergent. However, gauge 
invariance dictates a result proportional to F? v rather than At. This requires 2 powers of 
external momenta, so that the integrand must in fact converge like d 4 k/k 6 for large loop 
momenta k. In the limit mp e 3> m e ^ dm (relevant due to LEP and other collider/accelerator 
constraints), the contribution of momenta larger than m Fe is ~ l/m 2 Fe . The leading l/mp e 
term can thus be safely obtained by expanding the integrand in powers of l/mp e and 
keeping only the first term. 
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This corresponds to "pinching" the box with one F e into a triangle involving only 
electrons and an effective dm-dm-e-e coupling given bywith the real couplings a, b given 
by a + ib = c*c r . For this set-up, computing the cross-section is a loop-textbook exercise 
for which we find: 
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C is a function of m e and m dm given by the Passarino-Veltman scalar integral. For m dm > 
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m e , this function develops an imaginary part corresponding to the formation of a real 



e + e pair subsequently annihilating into 2 photons, and giving the largest contribution 
for masses above 1 MeV. 

For m dm <C m e , C behaves as [— I /(2m 2 ) + rn% m /(3mf)], so that both terms of the 
cross section behave as m dm / \m e rriF e ) 2 ■ This limit is relevant to estimate the effect of 
heavier particles than the electron in the loop. For example, the contribution of the 
r lepton could be significant if the corresponding couplings (a T ,b T ) are larger than ~ 
{m T /m dm ) x (a e ,b e ) x (m^/rnpj (with md m < m r ), i.e. if they scale at least like usual 
Yukawa couplings. Since an independent detailed analysis is required to check whether or 
not such couplings can pass particle physics constraints, we prefer giving a conservative 
estimate based on the electron contribution only. The latter cannot be turned off without 
losing the 511 keV line signal. It therefore constitutes a safe lower bound for assessing 
the detectability of the line at E y = m dm . 

Within the pinch approximation, the cross-section relevant for the origin of the 511 
keV emission is: 
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with f3 e = \J\ — rn 2 e /m dm) which indeed for b = reduces to the expression used^ for 
large m Fe . After careful comparison with SPI data, Ascasibar et al.^ found 



^ = 2.6 10-* (^f) 2 cmVs. 



The 77 annihilation cross-section is then also determined by this measurement in 
terms of the ratio of annihilation branching ratios: 
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As announced, this ratio cannot vanish, whatever the value of a/b, so that a minimum 
77 flux is guaranteed. As m dm approaches m e from above, the ratio increases like /?~ 3 for 
a pure scalar coupling (6 = 0) and like /J" 1 for an axial one (a = 0). The ratio decreases 
almost linearly with the dark mater mass for > 1 MeV. In the table below, we give 
typical values of the ratio rj for the most conservative case (i.e. a = 0, Z^ 1 ): 



m dm (MeV) : 


0.52 1 5 20 


r](a = 0) : 


8.8 10" 5 1.4 10" 5 3.610 _e 8.11Q- 7 



Notice that the simple guess^ applied to the case of decaying DM 
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increases instead of decreasing with m& m . For a typical mass of 10 MeV, this guess 
overestimates the monochromatic flux by a factor 635 with respect to our result (EqlTJ. 
As we will see in the next section, such a factor is crucial to the line observability. 
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Figure 1: Left: flux of the monochromatic E 1 = m dm line from a 8 degree cone around the galactic 
center. Right: significance of the monochromatic E y = m dm line above the continuum background for 
one year of observation with an ideal detector of 1 m 2 and a 10 -3 energy resolution. 



3 Detectability of the monochromatic line 

A few experiments have already scanned the energy range above the electron mass. The 
instruments on board of INTEGRAL for example have been designed to survey point-like 
objects as well as extended sources over an energy range between 15 keV-10 MeV. The 
instrument INTEGRAL/SPI itself is a spectrometer designed to monitor the 20 keV-8 
MeV range with excellent energy resolution. Therefore a legitimate question is whether 
or not the line E 1 = md m could have been (or could be) detected by the same instrument 
that has unveiled the 511 keV signal. This essentially depends on the ratio r\ as given 
above, and on the background. 

The 511 keV emission has been measured with a ~ 10% precision ^to be 



(J 511 ) = 6.62 x 10~ 3 ph cm 



2 s 1 sr 1 



inside a region that extends over 350° < / < 10° in longitude and |6| < 10° in latitude. 
If this emission originates from a NFW distribution of LDM species around the galactic 
center with a characteristic halo radius of 16.7 kpc, the signal from the inner 5° is found^ 
to be 

(/ 5 ii(5°)> = U 



x 10 2 ph cm 2 s 1 sr 1 



once the SPI response function is taken into account and the instrumental background is 
properly modeled. If the positron propagation is negligible, then the map of the 511 keV 
emission should correspond to that of the LDM annihilations. 

Within this approximation, we expect the spatial distributions of both the 511 keV 
and the two-gamma ray lines to be identical and their intensities to be related by the 
ratio 77: 

#511 (-^51l(^51l)) 



(.T 77 (0 77 )) = V 
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This expression is approximately valid as long as the angular radii # 77 and #511 of the 
regions monitored by the gamma-ray spectrometer are small. In what follows, the fraction 
/ of positrons forming positronium has been taken SI equal to 93%. This is in perfect 
agreement with the positronium fraction later derived^, i.e. fp s = 0.92 ± 0.09. The 
monochromatic line flux 

77 (< 6 I 77 ) = 7r0 77 (J 77 (6 I 77 )) 

has been plotted in Fig. [TJin the case of the LDM model with F exchange and assuming 
a NFW profile. The angular radius # 77 = 8° corresponds to the field of view of the 
satellite. For typical LDM masses in the MeV range, the expected flux is about three 
orders of magnitude below the claimed INTEGRAL/SPI line sensitivity^ ( w hich is 
about 2.5 x 10~ 5 ph cm" 2 s" 1 after 10 6 seconds). An unrealistic exposure of 30,000 years 
would thus be required in order to detect the E = md m line. When the LDM species 
is degenerate in mass with the electron, the flux is only a factor of 25 below the SPI 
detection limit (assuming a pure scalar coupling b=0). As long as the mass difference 
m dm — m e does not exceed 0.1 MeV, it is roughly comparable with the expected 478 keV 
line signal emitted by Novae^, that is about ~ 10~ 7 ph cm" 2 s -1 . 

SPI sensitivity is limited by the instrumental background that arises mostly from 
cosmic rays impinging on the apparatus and activating the BGO scintillator. On the 
contrary, the absolute sensitivity of an ideal instrument is purely limited by the gamma- 
ray continuum background. This emission has been recently estimated^ 



I BG (E) = 1.15 x 10" 2 TT L82 ph cm" 2 s" 1 sr" 1 



MeV" 



■ 



inside the central region that extends over 350° < I < 10° in longitude and \b\ < 10° in 
latitude. The energy E is expressed in units of MeV. 

We thus estimate the significance S = signal/ ^background for the LDM line to emerge 
above this background (assuming it is isotropic) to be 

7T0 




511 

with So the surface of the detector, T the exposure time, Ibg the above-mentionned 
continuum background intensity and AE the energy resolution. The significance S (dis- 
played as a function of the dark particle mass in Fig. Q] for a surface of 1 m 2 , an exposure 
duration of To = 1 year and an energy resolution of 0.1%) indicates that those values 
would theoretically allow to extract the minimal guaranteed signal computed at 3 stan- 
dard deviations above background for all relevant LDM masses below 30 MeV. 

There is nothing to be gained by narrowing the angular aperture # 77 because, for the 
assumed NFW profile, the signal increases linearly with this angular radius, as does the 
square root of an isotropic background. 

In contrast, note that the monochromatic line should be extremely narrow: its width 
is expected to be about a few eV which experimentally is very challenging if one compares 
it with the present SPI sensitivity that is about 10~ 3 at MeV energies. At lower energies, 
there are nevertheless instruments, e.g. X-ray CCD, bolometers, Bragg spectrometers 
which are able to resolve eV widths. A significant improvement on the resolution AE at 
higher energies would probably be necessary in order to reach a large enough significance 
and ensure detection. Indeed, an effective surface of lm 2 might be hard to attain in space. 
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Next generation instruments such as AGILE/ (super AGILE) or GLAST, which in 
principle could be more promising, will probably be limited by the energy range that they 
are able to investigate. Future instruments might nevertheless be able to see this line if 
their energy resolution and sensitivity are improved by a large factor with respect to SPI 
present characteristics. 

Maybe a better chance to detect this line is to do observations at a high latitude and 
a longitude slightly off the galactic centre. In this case, indeed, the ba ckgr ound should 
drop significantly (the density of dark clouds has been measured recently^) but the line 
flux may decrease by a smaller factor. 
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